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We derive a general formula of the reduced fidelity susceptibility when the reduced density matrix 
is 2 x 2 block-diagonal. By using this result and the continuous unitary transformations, we study 
finite-size scaling of the reduced fidelity susceptibility in the Lipkin-Meshkov-Glick Model. It is 
found that it can be used to characterize quantum phase transitions, implying that we can extract 
information of quantum phase transitions only from the fidelity of a subsystem, which is of practical 
meaning in experiments. 
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I. INTRODUCTION 

During the past few years, some important concepts 
in quantum information theory have been introduced to 
characterize quantum phase transitions (QPTs). For ex- 
ample, entanglement, which is one of the central concepts 
in quantum information theory, has been investigated ex- 
tensively in QPTs in various models, like Ising model 
[1-4] and Lipkin-Meshkov-Glick (LMG) model [5]. Re- 
cently, fidelity, which is another important quantum in- 
formation concept, has also been applied in characteriz- 
ing QPTs. The introducing of fidelity in QPTs is natural 
[6-23], since it's mathematically the overlap between two 
states, while QPTs are just dramatic changes in ground- 
state properties. However, the fidelity used in the study 
of QPTs depends computationally on an arbitrary yet 
finite small change of the driving parameter. To can- 
cel the arbitrariness, Zanardi et al. introduced the Rie- 
mannian metric tensor [15], while You et al. suggested 
the fidelity susceptibility [11]. The fidelity susceptibility 
then becomes an effective tool to study critical properties 
[15, 16] in many-body systems. 

It's noticed that all the above works are concentrated 
on the fidelity of the global ground states, and we 
may call this kind of fidelity susceptibility the global fi- 
delity susceptibility. However, in experiments, one al- 
ways probe the subsystem but not the whole system for 
practical convenience. Here we use the reduced fidelity 
[24] (also called partial fidelity in [25, 38]) susceptibility 
(RFS), which describes the fidelity susceptibility of a sub- 
system. In this work, first we derive a general formula 
of the reduced fidelity susceptibility when the reduced 
density matrix is 2 x 2 block-diagonal. Then, considering 
the LMG model, we show that the RFS can be used to 
characterize QPTs, and find that the scaling exponent is 
different from that of the global fidelity susceptibility. 

This paper is organized as follows. In Sec. II, we briefly 



review the concept of fidelity susceptibility, and give a 
general formula of RFS for a special but interesting case 
that the density matrix is 2 x 2 block-diagonal. Then 
in Sec. Ill, we introduce the LMG model [26]. in the 
isotropic case, we find that the critical behavior of RFS 
X in response to magnetic transverse field h as (h c — h)~ l 
in thermodynamic limit. While in the anisotropic case, 
by using the continuous unitary transformations (CUTs) 
[27-29], we find that the maximum of x over h diverged 
as iV 2 / 3 for an TV-spin system, and \h c — h\ in thermo- 
dynamic limit. Finally, we perform a numerical scaling 
analysis, and the results are well consistent with our the- 
oretical ones. 



II. REDUCED FIDELITY SUSCEPTIBILITY 

We first give a brief review on the concept of fidelity 
susceptibility. The Hamiltonian of a quantum system 
undergoing QPTs can be written as 

H(h) = Ho + hHj, (1) 

where Hi is supposed to be the driving term with control 
parameter h. The global fidelity is defined as F (h, S) = 
\{tpo (h) \<po (h + S))\, where \fo (h)) is the ground state 
of H (h), and 5 is a small quantity. The reduced fidelity 
is defined as the overlap between the reduced density 
matrix (RDM) p (h) of the ground state \ipo (h)). In the 
follows, we take p = p (h) and p = p(h + 5). Then the 
reduced fidelity is given by [30] 

F (h, S) = trvV/W/ 2 . (2) 

The corresponding fidelity susceptibility is defined as [7, 
11] 

r - 21nF ™ 



'Electronic address: xgwang9zimp.zju.edu.cn 



and then we could write F ~ 1 — x<5 2 /2. 
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In this papar we consider that the RDM is block- 
diagonal, 



P = 0ft, 



(4) 



i=i 



where ft 's are 2x2 semi-positive definite Hermitian ma- 
trices, since p is a density matrix. Now we introduce 
some useful formulas at first. Let A and B are arbitrary 
2x2 semi-positive definite matrices, then we have 



tr\ / A 1 / 2 BA 1 / 2 = y/tr (AB) + 2\/dct (AS), (5) 

and if A = B, it becomes 

tr(^ 2 ) = (tr^) 2 -2 dot A. (6) 

Take derivations of the above equation with respect to 
some variable h, we get 

tr (AA') = trAtiA' - d h (dot A) , (7) 

tr (AA") = trAtrA" - d 2 h (det A) + 2 det A', (8) 

where A' = d h A, A" = d 2 A and d h tr(A) =tr(A'). Now 
the fidelity can be written as 



1/2- 1/2 

QiQi 



= X v trft ^ + 2 V det ft ^ 



(9) 



i=l 



and recall that F ~ 1 — x<5 2 /2, the susceptibility x = 
S"=iXij with Xi corresponds to the 'susceptibility' of 
the i-th block in Eq. (4) . To obtain the susceptibility, we 
should expand the fidelity with respect to 5, and for ft 
ft + Qi W 5 + 5 2 g'l (h) /2 + (h 3 ) , we have 



tr (gg) ~ tr (g 2 ) + tr (ft/) $ + -tr (gg") , 

<5 2 

det g~detg + d h (det g) <5 + — d\ (det , 



(10) 



here we omit the subscript i for convenience. 

In the case that det g ^ 0, we have trg ^ since g is 
semi-positive definite. Then we get 



\J det (gg) ~ det g + -<9/j det g 
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9^ det g ■ 



(dh dot g)' 
2dctp 



(11) 



Take the above expression into Eq. (9) and with the help 
of Eqs. (6), (7) and (8) we obtain 



j ^2 

trvV /2 ft? 1/2 ^ trg + -tig' + -jtrg" 
+ ^Lactg>~(trg>f-^ dCt ^ 



8trg 



det (g) 



(12) 



If det g = but trg 7^ 0, we have det (gg) = 0. Moreover, 
since g is positive semi-definite, zero is the lower bound 
of dot g, which requires dh det g = and det g > 0. 
Thus we have 

tr(ft?) = (trg) 2 + trgtvg'6 
S 2 

+ y [trgtrg" - d 2 h (det g) + 2 dot £>'] , (13) 



and 



r r2 

trV 'g 1 / 2 gg 1 / 2 ~ trg + -trg' + — trg" 



8trg 



4detg' - (trg') - 2^ (detg) 



(14) 



In the last case that tig — 0, g is equivalent to a zero 
matrix, since g is Hermitian. Then tr(gg) = \J det (gg) = 
0, and F = 0. 

Conclude the above three cases, we get the 'suscepti- 
bility' for block Qi as 



1 I 1. m2 „ , , , det (ft 

; (trft) -4 det ft + 



4trft 
1 

4trft 




det (ft) 



(trftr-4det^ + 25 2 (detft 



r 



for trft 7^ 0, det ft 7^ 0, 

for trft 7^ 0, det ft = 0, 
for trft = 0, 



(15) 



where the terms of #trg'/2 and <5 2 trg"/4 in Eqs. (12) and due to tr(p) = 1, and ti(p') = tr (p") = 0. 

(14) are canceled in the final expression of the fidelity, _,. „ . , . , .... 

h mally, we consider a more special case that p is diag- 
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onal, the then susceptibility is obtained readily 

\2 



X 



(16) 



where Aj's are the nonzero diagonal terms. 



III. THE LMG MODEL AND ITS SCALING 
EXPONENTS OF RFS 

A. The LMG model and RFS 

The LMG model was introduced in nuclear physics to 
describe mutually interacting spin-1/2 particles, embed- 
ded in a transverse magnetic field. In the thermodynamic 
limit, it undergoes a QPT that is described by the mean 
field analysis [31]. Recently the finite-size scaling was 
studied by the 1/N expansion in the Holstein-Primakoff 
single boson representation [32] and by the CUTs [33, 34]. 
The Hamiltonian of the LMG model reads 



H =H + hHi 

= -A(l+ 7 ) (&-SI-N/2) 

-^(1-7) (S 2 + + S 2 _)-2hS z , 



(17) 



where S a = < T ia/2, with a a (a = x, y, z) the Pauli ma- 
trices, and S± = S x ±iS y . The prefactor 1/N ensures 
finite energy per spin in the thermodynamic limit. In the 
context, we set the parameters: A = 1, |7| < 1, h > 0. 
We take h > as the spectrum is invariant under the 
transformation h <-> —h. In addition, we only consider 
the maximum spin sector S = N/2 in which the lowest 
energy state lies. 

Now we consider a 2-body RDM of the LMG model 
[35] 



Pij 



fv+ u 

y y 

y y 

u w_ 



(18) 



in the standard basis 
a z \0) = -|0) and <r x \l) = 
elements reads 



{|00),|01),|10),|11)}, where 
|1), while the nonzero matrix 



N 2 - 2N + A(S 2 ' 



v± 



V : 



±4(S Z )(N-1) 



AN (N - 1) 
N 2 -A{S 2 ) _{S 2 X -Sl) 



AN (N - 1) 



N(N -I)' 1 



(19) 



where [A, B] + = AB + BA is the anti-commutator for 
operators A and B. The zero elements of result from 
the fact that the total spin and the parity are conserved 
quantities, i.e., 



[H,S 2 ] = 



N 



0. 



(20) 



It's noticed that is actually block-diagonal in the re- 
arranged basis { 1 00) , |11), |01), 1 10)}, and the two blocks 



are 



v + u 

U V- 



V V 

y y 



(21) 



With the help of Eq. (15), we can give the RFS explicitly 



2y 4(u++-y_) 
(v' + v_ + v + v'_ — 2u'i 



V '_) 2 +Au' 2 



+ 



(v+V- — u 2 ) 



(22) 



here we consider the case that det Q\ 7^ 0, and the follow- 
ing computations are based on the above formula. 



B. The isotropic case 

Firstly, we consider the isotropic case, 7 = 1, and the 
Hamiltonian reads 



H 



N 



(S 2 - S 2 - N/2) - 2hS z 



(23) 



which is diagonal in the standard eigenbasis {\S, M)} of 
S 2 and S z . For S = N/2 the eigenstates are 

E(M,h) = UM- f fY-? [ (l + h 2 ), (24) 



and the ground state is readily obtained when 

U (! : - i .. ". ..... ... . (25) 



f N/2 for h>\, 

''' ! = [N/2-R [N(l - h)/2] for < h< 1, 



where R(x) = round(a;). Then one can see level crossings 
exist at h = hj, where h j = 1 — (2j + 1) /N, between the 
two states \N/2, N/2-j) and \N/2, N/2-j - 1). In the 
thermodynamic limit, these critical points form a region 
of criticality. 

The elements of the RDM in ground state are readily 
obtained as 



v± = 



y = 



(N ± 2M ) {N - 2 ± 2M ) 
AN (N - 1) : 
(iV 2 - 4M ) 
AN(N-l) 



u = 0. 



(26) 



As N is very large, M (h < 1) ~ hN/2. With Eq. (22), 
we obtain the susceptibility in thermodynamic limit 



lim v I h > 1 — ^7 



2(l-/i 2 )' 



(27) 



Obviously, the asymptotic behavior of x as /i ^ 1 is 
1/ (1 — ft.). However, there is no QPT in its symmetric 
phase h > 1, because the ground state is independent of 
h. 
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FIG. 1: Fidelity susceptibility \ as a function of h with various 
system size N = 2 7 , 2 8 , 2 9 , 2 10 . The positions of their peaks 
approach to the critical point h c — 1 . 



C. The anisotropic case 

1. Spin expectation values 

Next we consider the anisotropic case, and the numer- 
ical results of the RFS as a function of h are shown in 
Fig. (1). We adopt the 1/N expansion method with 
CUTs that was used extensively by Dusuel and Vidal 
[33, 34], which corresponds to the large N limit. While 
the Holestein-Primakoff method is not suitable for our 
task since it could only give a first order correction in a 
1/N expansion. 



Here we firstly recall the CUTs introduced by Wegner 
[27] and independently by Glazek and Wilson [28, 29]. 
For a pedagogical introduction to this technique, one can 
sec [36]. The main idea of CUTs is to diagonalize the 
Hamiltonian in a continuous way starting from the orig- 
inal Hamiltonian H = H (I = 0). A flowing Hamiltonian 
is then defined by 



H(l) = C/t (l)H(Q)U(l). 



(28) 



where U (I) is unitary and I is a scaling parameter such 
that H (I — oo ) is diagonal. A derivation of the Eq. (28) 
with respect to / yields the flow equation 



d l H{l) = [r ] {l),H{l)], 



(29) 



where rj(l) = —U^diU is an anti-Hermitian generator. 
To obtain the expectation value of any operator on an 
eigenstate \ip) of H , one should follow the flow of the op- 
erator fi (I) = (I) H (0) U (I), by solving Eq. (29). For- 
tunately the results of the spin expectation values have 
been obtained by Dusuel and Vidal in [33, 34], and here 
we'll compute the scaling behavior of the derivatives of 
these values. 

Firstly, we consider the system size N is very large, 
and the matrix elements are rewritten as 



"±. 

y 



N 



1 i (ff) ■ (S*) 
4 N 2 

4 N 2 



(s 2 x ) - (S 2 y ) 



N 2 



(30) 



The spin expectation values appeared in the above ex- 
pressions can be solved by the CUTs with 1/N expan- 
sion. For symmetry phase (h > 1), we have 



2(3,) 
N 

MS 2 X ) 

N 2 

MS 2 ) 

N 2 

MS 2 ) 

N 2 




(31) 



where G = G (h, 7) = (h — 1) (h — 7). Here we do not 
present P ? (l) = P ? (l) (ft, 7) and Q® = (h, 7) (i = 
1,2,3 and £ — z,xx,yy,zz), which are polynomials of 
h and 7, whereas of little meaning for computing the 
scaling exponents. For more details, you can refer to 



the appendix part of [34]. It's noticed that, the above 
expressions can be written in the form 

1> N (h, 7) = $^ (h, 7) + 4>f S (h, 7) , (32) 
where the superscripts 'reg' and 'sing' stand for regular 
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and singular respectively. A nonsingular contribution is 
understood to be a function of h which is nonsingular at 
h = 1, as well as all its derivatives. Take 2(S Z )/N for 
example, the regular part is 1 + l/N and the remain- 
ing forms the singular part. As h approaches to 1, the 

(i) 

terms involving s are small compared to the terms 
involving Pg^'s by a factor G (h, 7), hence we could only 

(i) 

consider the terms involving P£ s. 



2. Finite-size scaling 

Here we show how to derive the finite-size scaling expo- 
nents of the spin expectation values and their derivatives, 
and take 2(S Z )/N for example, 




log 2 /V 



2(S Z ) 1 1 
x 1 = H 1 

N NG 1 / 2 



7) 2 Pi 3) 



N 



p(i) 



(1 



(2) 



NG 3 / 2 



(1 



(NG 3 / 2 



O 



[NG 3 / 2 



(33) 



where the singular part (terms after 1 + l/N) can be 
written in the form 



2(S Z ) 
N 



1 



NG(h,i) 



1/2 



NG(h n f 2 yl 



(34) 

where ($ = S z , S 2 , S 2 , S 2 ) is a scaling function 
for these spin expectation values. While in fact that 
there can be no singularity in any physical quantity in 
a finite-size system, and the critical point h c = 1 only 
for thermodynamic limit N — > 00. This implies that 
the singularity of G (h, 7) ' has to be canceled by 



the one of J^V. 



NG(h,j) 3/2 ,j 



Thus one must have 



^"s z ( x :l) ~ x 1 ^ 3 5 which in turn implies the following 
finite size scaling: 



FIG. 2: Maximum susceptibility \m as a function of system 
size N. We can see that the numerical results approach to 
the solid line with slope 2/3 as the system size increases. 



x 10" 



I 




0.05 

rf 3 (h-h) 



2{S Z 



N 



h=l 



N 2 / 3 ' 



(35) 



Immediately, one can obtain the asymptotic form of all 
the spin expectation values 



2(5,) 



N 

MS 2 ) 



1 



N 2 

Ms 2 y ) 



N 2 
MS 2 ) 



N 2 



h=l 



h=l 



h=l 



h=l 



1 

TV 



N 2 / 3 '' 



a (0) 
N 2 / 3 ' 

a (0) 
u yy 

AT4/3' 

2 

1+ N 



a (0) 
N 2 / 3 ' 



(36) 



where a^°''s (^ = z,xx,yy,zz) are all constants depend- 
ing on 7. Then take the first-order derivatives of Eq. (31) 



FIG. 3: Finite size scaling is performed. The susceptibility 
X is considered as a function of the system size TV and the 
parameter h, take the form iV" (h — h m )- Here the exponent 
v = 2/3 is determined analytically. It's noticed that the data 
dose not collapse at one line exactly since the system sizes are 
not large enough. 



with h, one could find similar scaling functions with 
Eq. (34). Here we also take 2(S Z )/N for example, 



( 9 2(S Z ) 
\dh N 



NG (h, 7) 



3/2 ' 



NG(h,~/) 3/2 n 



(37) 

where ^$ is a scaling function for the derivatives of spin 
expectation values, and then we find the finite size scaling 



d 2{S Z ) 



dh N 



(38) 



h=l 
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The scaling form of the other derivatives are 



dh N 2 
d MS 2 ) 



dh N 2 
d 4{S 2 ) 



dh N 2 



h=l 



h=l 



h=l 



a (1) 

a (1) 
a yy 

AT2/3' 
«(!) 



(39) 



where a^'s (£ = z,xx,yy, zz) are constants depending 
on 7. As we can see that, except for 4(Sy)/N 2 , the other 
first-order derivatives are all independent of N. Then 
with the help of Eq. (22), we find that the maximum 
RFS xm =x(hm,N,i) is 



Xr 



(°») S 



N 



(40) 



for large A/', and here we just present the divergent term. 
It's noticed that a^} should be less than —2 to ensure 
the matrix element y > 0, thus Xm > 0. Then we have 



Am In N + const. 



(41) 



where the constant only depends on 7 and the scaling 
exponent A n approaches to 2/3 as N increases, which 
is verified numerically, and An — 2/3 in thermodynamic 
limit. The numerical comparisons are shown in Fig. (2). 
While in the broken symmetric phase (0 < h < 1), we 
can derive the same scaling exponents [34] . However, for 
global fidelity susceptibility, the scaling exponent is 9/7 
[19]. 

Then if we cancel N in Eq. (37), with similar steps, 
we can get the relation between the susceptibility x an d 
77 = h — h c in thermodynamic limit, 



lnx (h, 7) = Ah In \h — h c \ + const. 



(42) 



where Ah approaches to —1 as h goes to h c , and the 
constant depends on 7. Therefore we could take the form 
of the susceptibility for finite size as 



x(h,N,j) = 



N- 2 /z + B{h~h m y 



(43) 



To study the critical behavior around the phase 
transition point, we could perform the finite scal- 
ing analysis. According to the scaling ansatz [37], 
the susceptibility is a function of N u (h — h m ). In 
the case of logarithmic divergence, it behaves as 
X(h m ,N)/ X (h,N) ~ Q[N v (h-h m )], where the 
function Q(x) « In a; for large x is universal and does 
not depend on system size N. Hence with Eqs. (41) 
and (43), we determine the exponent v = 2/3, which is 
confirmed numerically, as shown in Fig. (3). However, 
the curves for different system sizes does not collapse to 
a single one exactly, since the system sizes are not large 
enough. 

IV. CONCLUSION 



In summary, we have investigated the RFS in the sec- 
ond order quantum phase transition of the LMG model. 
For the case that p is block-diagonal in 2 x 2 matrices, we 
derive a general formula for RFS. Then with the CUTs 
and the scaling ansatz, the critical exponents, includ- 
ing the finite-size scaling exponents of the RFS are ob- 
tained analytically, and confirmed numerically. Our re- 
sults show that, the RFS undergoes singularity around 
the critical point, indicating that the RFS can be used 
to characterize the QPTs. And it's suggested that we can 
extract information of the QPTs only from the fidelity of 
a subsystem, without probing the global system, which 
is of practical significance in experiments. It is also inter- 
esting to study finite-size scaling of RFS in other models 
such as quantum Ising model, which is under considera- 
tion. 
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